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Abstract 

There are at least three models of discrete-time quantum walks (QWs) on graphs 
currently under active development. In this work we focus on the equivalence of two 
of them, known as Szegedy’s and staggered QWs. We give a formal definition of the 
staggered model and discuss generalized versions for searching marked vertices. Using 
this formal definition, we prove that any instance of Szegedy’s model is equivalent to an 
instance of the staggered model. On the other hand, we show that there are instances 
of the staggered model that cannot be cast into Szegedy’s framework. Our analysis also 
works when there are marked vertices. We show that Szegedy’s spatial search algorithms 
can be converted into search algorithms in staggered QWs. We take advantage of the 
similarity of those models to define the quantum hitting time in the staggered model 
and to describe a method to calculate the eigenvalues and eigenvectors of the evolution 
operator of staggered QWs. 


1 Introduction 

Quantum walks (QWs) are being used with many flavors with the goal of understanding 
quantum systems and building quantum algorithms for quantum computers. Research in 
QWs started with the coined discrete-time QW [T] proposed with the goal of displaying a 
quantum system with features strikingly different from classical systems. From a different 
viewpoint, Meyer [2] proposed a non-trivial example of a quantum cellular automata which 
can be somehow considered the starting point of staggered QWs. By quantizing continuous 
classical Markov chains, Farhi and Gutmann [3] proposed a continuous-time version of QWs, 
with clear implications to the area of quantum computing. Inspired by classical discrete 
Markov chains, Szegedy [3] proposed a coinless discrete-time QW model and was able to 
provide us with a natural definition of quantum hitting time. 

‘Corresponding author: portugal@lncc.br 


1 



Szegedy’s QWs are obtained by quantizing classical discrete Markov chains described by 
some transition matrix. Szegedy [1] also developed QW-based search algorithms inspired 
by a previous known coined-based search algorithm [5]. On ergodic Markov chains, it is 
possible to detect the presence of a marked vertex at a hitting time that is quadratically 
smaller than the classical average hitting time |^. Szegedy’s model was also used for 
the searching problem laEj, which aims to find the location of a marked vertex, and for 
searching triangles [H]. Some references highlighted the power of the method, for instance 
see Ref. [9]. The fact that Szegedy’s model is a QW on the edges of a bipartite graph 
rather than on the vertices is emphasized in mM- In this paper we investigate this idea 
by employing the line graph of the bipartite graph. 

Meyer [2] circumvented the triviality of quantum cellular automata by relaxing the ho¬ 
mogeneity condition. He proposed an one-dimensional quantum cellular automata driven by 
staggered unitary operators, which can be transposed as a coinless QW on the line or cycle. 
The connection of this model with the staggered lattice fermion formalism was analyzed 
in mi for the two-dimensional case and applied for spatial searching in [12]. The staggered 
model was rediscovered by Falk |I3|, who suggested a simple method of obtaining the evo¬ 
lution operator by splitting the vertices of the graph into disjoint polygons that tessellate 
the two-dimensional lattice. Falk also proposed a searching method, which was used by m 
to prove analytically that the staggered model hnds a marked vertex in time 0{yJN log Ai) 
for a two-dimensional lattice with N vertices matching the performance achieved by other 
forms of QWs on this graph. The detailed dynamics of the one-dimensional staggered 
model was analyzed in m and the moments of the probability distribution function in m- 
The connection between the coined and the staggered one-dimensional QWs was analyzed 
in nziiiaEii. Further analysis on this connection is an important research topic. 

In this paper we formally define the staggered QW model on graphs by using tessel¬ 
lations so that each tessellation covers all vertices with non-overlapping polygons and the 
tessellation union covers all edges. The maximal cliques play an essential role in building 
the tessellation because adjacent vertices are reachable for the walker’s hopping. If two or 
more maximal cliques share a common vertex, each clique must be inside a different tessel¬ 
lation. The evolution operator is a product of orthogonal reflections that have a one-to-one 
correspondence with the tessellations producing a diffusion through the maximal cliques in 
a staggered way. After the action of the evolution operator with a localized initial condi¬ 
tion, the wave function spreads to all vertices reachable by the graph structure, which is 
the clique with the starting location and the adjacent maximal cliques for the case with two 
reflections. We also discuss a generalized form of staggered QW using partial tessellations 
which can be used for spatial search algorithms in the sense that the walker will search 
for the vertices inside the missing polygons. The complete dynamics is obtained by apply¬ 
ing the evolution operator recursively starting with some initial state and by performing a 
measurement at the end. 

After formally defining the staggered model, we show that Szegedy’s QWs on a bipartite 
graph F are staggered QWs on the line graph of F with a restricted form of tessellation. 
The restriction is: The polygons must share only one vertex. The converse is also true. 
Any staggered QW with tessellations that share one vertex in the polygon intersections can 
be cast into an extended version of Szegedy’s model with complex amplitudes. This result 
sheds light on the structure of Szegedy’s model and explains the meaning of the fact that a 
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Szegedy’s QW is a hopping on the edges of the bipartite graph rather than on the vertices. 
We also show that there are generalized staggered QWs that reproduce the searching model 
proposed by Szegedy. In this case we need to use partial tessellations and the walker looks 
for vertices inside missing polygons, which are the ones returned by the measurement. On 
the other hand, if we allow two or more vertices in at least one polygon intersection, we 
define instances of staggered QWs that cannot be cast into Szegedy’s framework. We give 
special attention to those instances. 

Falk |13] proposed an evolution operator for a searching model on the two-dimensional 
lattice by interlacing a reflection around the marked vertices with the orthogonal reflections 
generated by the tessellations. We use this searching model to define the hitting time on 
staggered QWs on finite graphs by generalizing Szegedy’s definition and we describe how 
to obtain the eigenvalues and eigenvectors of the evolution operator using the discriminant 
matrix generated by the inner product of the polygons of the tessellations. We show that 
the evolution operator with marked vertices can be written as a product of two orthogonal 
reflections, and if an entire polygon is marked, the walker is not able to find it. 

The structure of the paper is as follows. In Sec. we define the staggered QW model 
and discuss generalizations that are useful for searching algorithms. In Sec. we review 
Szegedy’s QW model with the goal of finding the connection with the staggered QW model. 
In Sec. 1^ we show that Szegedy’s model is a restricted form of staggered QWs and we 
characterize when staggered walks cannot be cast into Szegedy’s framework. In Sec. [^we 
define the hitting time on the staggered model using reflection operators around the marked 
vertices. In Sec. we analyze the eigenvalues and eigenvectors of the evolution operator 
for spatial search algorithms based on reflections around the marked vertices. In Sec. we 
draw our conclusions. 


2 Staggered Quantum Walks 

Inspired by Falk’s paper m, we define the staggered quantum walk on a graph F with 
N vertice^by using two (or more) independent graph tessellations. Each tessellation uses 
non-overlapping (non-planar) polygons of adjacent vertices which need not to have the same 
shape and a polygon may contain only one vertex. The set of polygons of each tessellation 
need to cover all vertices of the graph and the tessellation union must cover all edges. Two 
polygons in different tessellations necessarily overlap and some (or all) intersections may 
contain more than one vertex. Each polygon defines a unit vector in the Hilbert space 
by superposing the vertices in the polygon with nonzero amplitudes. Examples of 
tessellations are depicted in Eigs. [^and[^ 

The staggered quantum walk with two tessellations called a and f3 is defined by the 
evolution operator 

U = UiUo, (1) 

where 

m—1 

Uq = 2'^^\ak){ak\ — I, ( 2 ) 

k=0 

can be infinite. 
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( 3 ) 


n—1 

Ui = 2Y,\Pk){Pk\-I, 

k=0 

and m and n are the number of polygons in each tessellation, and 


|aft) 


(4) 

|A> 

^ ^ bk,k' 1^)5 

k'e/Sk 

(5) 


where ak,k' and bk,k' are nonzero complex amplitudes of the unit vectors la^) and \Pk) 
in T-L^. Operators Uq and Ui are unitary and Hermitian [Uqi = l) because each set of 
polygons is non-overlapping {{ak\ak’) = (/3fc|/3fc') = 6kk')- It is straightforward to show 
that 

U\k) = 4 Dl,,„al,k\Pk") -2Y,bhk\Pk") - 2 + \k), (6) 

k\k" k" k' 

where Dyk" = {o^k'\Pk")^ k' runs from 0 to m — 1, k" from 0 to n — 1, and k from 0 to 

At- 1. 



(a) Cycle with 2-site poly¬ 
gons 



(b) Glued triangle tree 


Figure 1: Examples of tessellations for the 10-cycle and a glued triangle tree with 22 vertices: 
(a) Depicts tessellations with polygons of size 2 for the cycle and (b) depicts tessellations with 
polygons of size 3 (except for two vertices) for the glued triangle tree. Ug is associated to the blue 
tessellation (solid line) and Ui is associated to the red tessellation (dashed line). 

In order to give an alternative form of the definition of staggered QWs, let us define the 
notion of orthogonal reflection of a graph. 

Definition 2.1 A unitary and Hermitian operator U is called an orthogonal reflection of 
a graph if the eigenvectors of U associated with eigenvalue 1 have non-overlapping nonzero 
entries and the sum of those eigenvectors has no zero entries in the orthonormal basis 
associated to the vertices of the graph. 


4 






















This definition is basis independent because a basis change keeps the orthogonal reflection 
invariant and changes the orthonormal basis associated to the graph. A staggered QW on 
a graph T can be alternatively defined by a unitary evolution operator that is a product of 
orthogonal reflections of T that covers the edges of the graph. Given an orthogonal reflection 
in the basis of associated with the vertices of T, it is possible to build a disconnected 
A'-graph as a disjoint union of cliques, where each clique is associated to an eigenvector 
with eigenvalue 1 and the clique is formed with the vertices that have nonzero amplitudes. 
Each clique defines a polygon and the set of polygons is the tessellation associated with this 
orthogonal reflection. The product of orthogonal reflections defines the final graph after 
the union (and merging) of the edges of all cliques. This final graph must be equal to T. 
The union of the tessellations does not have disposable edges in the sense that every edge 
is inside a polygon of some tessellation. 

It is tempting to remove the requirement that the polygons must be cliques. The removal 
of this requirement is troublesome because after one application of the unitary operator Uq 
(or Ui) with non-clique polygons, the walker can go to non-adjacent vertices. In quantum 
walk models, we expect that the walker only hops to adjacent vertices, and successive 
applications of shift operators move the walker to distant locations. Our definition does 
not permit the tessellations used by Falk in Ref. m, which were also used by Ambainis et 
al. in Ref. mi and Patel et al. in Refs. mm- Those papers go against the requirement 
described in Ref. which states that a general form of quantum walks must respect the 
structure of the graph. The tessellations used in those references in fact refer to the graph 
depicted in Fig.[^ For the two-dimensional lattice (with degree 4), it is necessary to use at 
least four tessellations because every vertex belongs to four maximum cliques of size two. 
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Figure 2: Example of tessellations with polygons of size 4 on a non-planar regular graph of degree 6. 
Uq is associated to the blue tessellation (solid line) and Ui is associated to the red tessellation 
(dashed line). 

It is also tempting to remove the requirement that all edges must be inside of a polygon. 
The removal of this requirement is troublesome because it generates a logical indetermi¬ 
nacy, similar to one in the discussion about the use of non-clique polygons in the previous 
paragraph. If an edge is not inside a polygon, it means that it can be removed generating 
no change in the dynamics. In this case non-isomorphic graphs would have exactly the same 
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dynamics and would lead to confusion and logical problems. 


2.1 Generalized Staggered Quantum Walks 

Usually quantum walk models are defined with no marked vertices. In this case the dynamics 
has no special vertex and the amplitude associated with some specific vertex will not increase 
above average unless the initial condition is a special one. For instance, the coined model 
defined in Ref. |19) has those characteristics. The standard definition must be changed in 
order to mark a vertex. In the coined case, the coins for the marked vertices are different 
from the coins for non-marked vertices pO] . 

We can extend the staggered quantum walk definition by removing the requirement that 
each tessellation must cover the entire graph. It is allowed to use partial tessellations so 
that the partial-tessellation union covers all vertices (not necessarily all edges). The vertices 
that are not inside all tessellations are the marked ones. An equivalent way of obtaining the 
same extension is by allowing zero amplitudes in the entries of polygons Ok and I3k ■ The goal 
of this generalization is to define search algorithms. For example, Grover’s algorithm [21] 
can be seen as a generalized staggered quantum walk on the complete graph in the following 
way: Tessellation a has polygons of size one for each non-marked vertex and tessellation 
(5 has only one polygon with all vertices. The vertices that do not belong to tessellation a 
are the marked ones. The evolution operator of this generalized staggered model is equal 
to the one used in Grover’s algorithm. 


3 Szegedy’s Quantum Walks 


In this section we review Szegedy’s QW model with the goal of hnding the connection with 
the staggered QW model. Gonsider a connected bipartite graph F(X, Y, E), where X, Y are 
disjoint sets of vertices and E is the set of non-directed edges. Let 


0 A 
0 


(7) 


be the biadjacency matrix of F(X, Y, E). Using A, define P as a probabilistic map from X 
to Y. Using A'^, define Q as a probabilistic map from Y to X. If P is an m x n matrix, Q 
will be an n X m matrix, both are right-stochastic with the property that each row sums to 

1 . 

To define a quantum walk on the bipartite graph, we associate it with a Hilbert space 
'}{mn _ ^ where m = |X| and n = |y|, the computational basis of which is 

||x,y) : X G X, y G y|. Szegedy’s quantum walk is defined by the evolution operator 


where 


w = Ri Ro, 

(8) 

Ro = \4>x){4>x\ - I, 

X 

(9) 

Pi = 2'^\^Py){^jJy\- I, 

(10) 


y 
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and 


\(t>x) 

y&Y 

(11) 


xex 

(12) 

It is straightforward to show that 



W\x,y) = 

y'& 

'x\'lpy'} ~ ‘^\/Qyx\'lpy) ~ ‘^^/Pxy\<Px} + 2/); 

(13) 

where Cyx = For a detailed review see [22]. 



4 Equivalence between Szegedy’s and Staggered QWs 

In this section we show that any instance of Szegedy’s model is equivalent to a staggered 
quantum walk version. The converse is not true. We show that some instances of the 
staggered quantum walk model are equivalent to Szegedy’s model and we characterize when 
staggered walks cannot be cast into Szegedy’s framework. 


4.1 Szegedy’s Walks are Staggered Walks! 


Although Szegedy’s model is dehned on the Hilbert space associated with a bipartite 
graph T{X,Y, E), the dynamics takes place in the subspace spanned by the edges of T. 
W acts trivially on vectors \x,y^ that do not belong to the bipartite graph and the initial 
condition does not include those edges. In this section we show how to define a staggered 
QW based on generic stochastic matrices P and Q equivalent to Szegedy’s QW. 

Let N be the number of edges of T. Define two sets of polygons in the line graph L(r) 
by 


\ax) = ^V^/(x,y)), 

(14) 

y&Y 


\l^y) = X] 

xex 

(15) 


where / is a bijection between the edge set Li of T and the labels k of L(r), which has 
associated the Hilbert space with the computational basis { |^), = 0, ...,N — l|. Each 

polygon Ux is a clique and the union of polygons ax tessellates L(r). The same is true for 
polygons Py. 

Suppose that f{x, y) = k, where {x, y} G E. Then 


( n—1 \ / m—1 \ 

~ - Ij \fix,y)). (16) 

y'=0 j V x'=0 j 

Using that {ax'\fix,y)) = y/pPy5xx' and {ay’\f{x,y)) = ^/^Syy>, we obtain 

f/| ^ ^ l/^y') ~ ‘^^/^yx\l^y') ~ 2 | ) P \ f {x, y)') ■ (H^) 
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Eqs. (© and ( [I^ are essentially the same results if we apply bijection / to the labels (x, y) 
of the kets |x,yy^ of Eq. (13). This proves 


Proposition 4.1 Any instance of Szegedy’s quantum walk model on a bipartite graph T is 
equivalent to a staggered quantum walk version on the line graph ofT. 


4.2 Which Staggered QWs are instances of Szegedy’s Model? 

Consider a staggered quantum walk on a Wgraph E' with the following restriction: The 
intersections of polygons belonging to different tessellations contain one vertex. If polygons 
Ofc and Pk' share one vertex, the ket label of this vertex in Szegedy’s model will be \k, k'^. 
Because all vertices must be in both tessellations, this labeling method establishes a bijection 
f between the computational basis of used in the staggered model and the set of labels 
of the form \k, /c'), that will be used to build an instance of Szegedy’s model on a bipartite 
graph T, the line graph of which is r'0 By construction, the set of labels of the form 
\k,k''^ belonging to polygon shares the same value of k with different values of k'. After 
applying the bijection f to the kets of polygon we obtain a vector \oik) in (m, n 

are the number of polygons in tessellations a, (5) given by 

|dfc) = (18) 

k' 

with the same amplitudes of the original polygon ak, where k' runs over the subindices of 
polygons fdk' that have overlap with that is, fdk' H a*, / 0. Analogously, after applying 
the bijection /' to the kets of polygon Pk', we obtain a vector \Pk') in 4 ^™"' which has same 
value of k' as in the second slot of the kets, and it is given by 

1^*:') = '^bk',k\k,k'), (19) 

k 

where k runs over the subindices of polygons ak that have overlap with Pk'-, that is, 7 ^ 

0 . 

We must consider a second restriction regarding the amplitudes and bk\k of the 
staggered model. We assume that ak^k' and bk',k are non-negative real numbers obeying the 
constraints Ylk'^kk' ~ ^ Ylk^k'k ~ ^hat is, \ak) and \Pk') are real unit vectors. 
After imposing this restriction, such instances of the staggered quantum walk model can 
be put into the standard Szegedy’s framework. In fact, the evolution driven by U, which is 
defined using \ak) and \Pk) with the restricted version of the amplitudes, is equivalent to 
the evolution driven by W using \ak) and \Pk’) instead of \4>x) and The equivalence 

is brought out by applying bijection f to the labels of the kets of U\k^ given by Eq. ^ 
and comparing with IT|/'(fc)). The second restriction (real amplitudes) can be eliminated 
by a straightforward extension of Szegedy’s model, which is addressed in the next section. 

The instances of the staggered quantum walk models that cannot be cast into Szegedy’s 
framework with the above method are those that have at least one overlapping polygons 
of tessellations a and P with more than one vertex. To prove this fact, let |ji) and |j 2 ) 
be kets associated with two vertices in the overlap between polygons Ok and Pk'- There 

^Given a line graph F', there is only one bipartite graph F such that I/(F) = F' | 23| . 
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is a contradiction because |ji) and |j 2 ) must be mapped to kets of the form \k,k'^ both 
sharing the same label k and k'. In fact, in order to build bijection f, we need to map 
|ji) to |A:, G and |j 2 ) to |A:,G sharing the same k because both belong 
to polygon Ofc. We have to take k'l / k '2 because different vertices must have different 
labels. Besides, |ji) and \j 2 ) belong also to polygon Pk', then labels k[ and must be the 
same, which is a contradiction. Therefore, there is no bijection between the computational 
basis of T-L^ and the set of kets of the form \k,k'^ with the property of sharing the same 
k for polygons of tessellation a and simultaneously sharing the same k' for polygons of 
tessellation jS. 

The simplest non-trivial examples of staggered QWs that cannot be cast into Szegedy’s 
framework are obtained on the graph with only two connected vertices {N = 2). Take 
tessellations |a) = (|o) + |l))/-\/2 and |/3) = jV'), where 

9 9 

1^) = cos - |o) + sin- |l), 

for some angle 9 that is not a rational multiple of vr. The evolution operator is 


U = 


sin 6 cos 9 
— cos 9 sin 9 


which is non-trivial because there is no integer r such that 17’' = I. On the other hand, any 
Szegedy’s quantum walk W on any graph with two edges is trivial because = I. Notice 
that we do not need to go back to Szegedy’s framework to prove this result because we can 
use the equivalent staggered versions. This proves 


Proposition 4.2 There are instances of the staggered quantum walk model with two tessel¬ 
lations that eannot be cast into Szegedy’s quantum walk framework. Those instances have 
at least one polygon intersection with at least two vertices in common. 


For example, the glued triangle tree of Fig. [^cannot be cast into Szegedy’s framework 
because it is not a line graph. Even if it was the line graph of a bipartite graph, it would not 
be possible to obtain an equivalent Szegedy’s QW because there are polygons intersections 
with two vertices in common. 

An alternative proof of the above proposition in terms of graph theory is as follows. 
If the quantum walk graph F is not a line graph, by inspection we verify that the nine 
forbidden Beineke induced subgraphs [2Tj require either tessellations with two vertices in 
common or at least three-tessellations. If graph F is a line graph, then it has a Krausz 
partition [25]. If the Krausz partition is two-colorable, F is a line graph of a bipartite graph 
and can be reduced to Szegedy’s framework by using Eqs. (18) and (19). If the Krausz 
partition is not two-collorable, F requires either tessellations with two vertices in common 
or at least three-tessellations. 


4.3 Extension of Szegedy’s Model 

Szegedy’s model can be easily extended by allowing phases in vectors 

|x,2/), (20) 

y& 
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( 21 ) 


\i^y) = X] 

x&X 


Entries are non-negative real numbers obeying the constraints ~ 

as in the original model. It is straightforward to show that W is unitary in the extended 
Using the same algorithm of Sec. 4.2, we conclude that any staggered quantum walk 


case. 


with the property that the intersections of polygons belonging to different tessellations 
contain one vertex is equivalent to an instance of the extended version of Szegedy’s model. 
This proves 


Proposition 4.3 Any instance of the staggered quantum walk model with two tessellations 
which has only one vertex in common in each polygon intersection can he east into the 
extended version of Szegedy’s quantum walk framework. 


Notice that a staggered quantum walk model on a graph T with two tessellations which 
have only one vertex in common in each polygon intersection has a two-colorable Krauz 
partition. Therefore T is a line graph of a bipartite graph. 


4.4 Searching in Szegedy’s Model 

Szegedy’s model can be used to search for a vertex in a set of marked vertices in some 
graph T. The strategy is an extension of the classical method, which uses random walks 
on the same graph. The key concept in the classical case is the hitting time, which is the 
average time to hit a marked vertex for the first time using a random walk with stochastic 
(or transition) matrix P after specifying some initial condition. There is a clever way to 
find the classical hitting time using the following idea: Convert the original graph T into a 
new directed graph T' (with a new stochastic matrix P') by removing the edges that leave 
the marked vertices. Marked vertices are converted into sinks. The hitting time obtained 
using P' is the same using P because as soon as the walker hits a marked vertex using 
some edge that comes from a non-marked vertex, the walker needs not to go ahead. The 
advantage of using P' is that it does not have eigenvectors associated with eigenvalue 1 (the 
spectral norm is smaller than 1), that is, / — P' is an invertible matrix. We need to invert 
I — P' to find the hitting time if we want to avoid the calculation of the limiting probability 
distribution. 

Szegedy proposed a quantum version of this procedure. Let r{X,E) be the original 
graph, where X is the set of vertices and E the set of edges. Define r(X, X', E') as a bipartite 
graph obtained from r(X, P) by duplicating X and by converting edges {xi,Xj} G E into 
{xi,x'j} G E' . Until this point, no vertex has been marked and the quantum walk described 
in Sec. can be used taking P = Q. As before, define r'{X,X',E') as a directed bipartite 
graph by removing the edges of r(A, X' ,E') that leaves the marked vertices of X and X'. 
Add new non-directed edges connecting a marked x with its corresponding copy x' for all 
marked vertices. The quantum walk described in Sec. [^can be used taking P' = Q' , where 
P' and Q' are the new stochastic matrices of r'(A, X', E'). With this framework and taking 




V xy 


( 22 ) 
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as initial condition in Szegedy showed that the detection problem on r'(X, X' ,E') 

can be solved with a quadratic speedup compared to the time complexity of the same 
problem using random walks with symmetric and ergodic stochastic matrix P in r(X, E)^. 

The method of converting Szegedy’s QWs into staggered QWs described in Sec. |4.1| 
can be used in the searching context after some minor modifications. Let L(r) be the line 
graph of T{X^X\E') and define polygons Ux and j3y as the ones in Eqs. (14) and (15) for 


all X and y that do not belong to the set of marked vertices. Those sets of polygons define 
partial tessellations of T(r) which is used to define a generalized staggered QW equivalent 
to Szegedy’s QW on T'{X,X'^E'). The initial condition in l-L^ must be 


l^o) = 


(23) 


xy 


where / is the bijection from the E to the labels of T(r). 

The last ingredient in Szegedy’s QW is the measurement, which is performed in the 
computational basis of X. At this moment, the copy X' is traced out. If the measurement 
is performed at the right time, which is the quantum hitting time, the probability of finding 
a marked vertex will be high enough. The version using generalized staggered QWs is 
obtained as follows. For x G A define projectors 


ria; = ^ \x'){x' 
x'eax 


(24) 


Notice that Yhx = I- By measuring observable 


O = ^ xn,, (25) 

x&X 

we obtain some value xq G X, which is expected to be one of the marked vertices if the 
measurement is performed at the hitting time. In the staggered version, we are looking for 
the missing polygons of the (partial) a-tessellation. This proves 


Proposition 4.4 Szegedy’s searching framework on a bipartite graph T can he simulated 
by instances of generalized staggered quantum walks on the line graph o/T. 


5 Hitting Time on Staggered Walks 

We have showed that any instance of the extended Szegedy’s model can be cast into an 
equivalent staggered quantum walk. The converse statement is not true. There are in¬ 
stances of staggered walks that are not equivalent to any instance of the extended version 
of Szegedy’s model. We take advantage of the similarity of both models to define a notion 
of hitting time on the staggered model, which is specially useful on those staggered walks 
that cannot be cast into Szegedy’s framework. 

Let M be the set of marked vertices. Inspired again by Falk’s paper m, the evolution 
operator Um for searching a marked vertex is defined by 

Um = RmUiRmUo, (26) 
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where 


( 27 ) 


Rm = 2 — I. 

m^M 

Rm is a orthogonal reflection and flips non-marked vertices. 

Definition 5.1 The quantum hitting time Hm of a staggered quantum walk starting from 
the initial condition |V’(0)) i-s the smallest number of steps T such that 

^EI||«*)>-|</>(0 )>|t>l-M, (28) 

t=0 

where \ilj{t)) = (C/m)*|V'(0)) ond N is the number of vertices. 

The hitting time depends on the initial condition. An important one in this context is 
the uniform distribution because the value 1 — |M|/Ai is the distance between the uniform 
distribution over all vertices and the uniform distribution over the marked vertices. 


6 Searching using Staggered Walks 

Rm can be interpreted as an operator associated with a partial tessellation with \M\ poly¬ 
gons each one with one marked vertex. There is a simpler way to analyze the evolution 
operator (26) by using the fact that RmUiRm is an orthogonal reflection. 

Proposition 6.1 The evolution operator Um for searching a marked vertex using a stag¬ 
gered quantum walk is a product of two orthogonal reflections. 

Proof Define U[ = RmUiRm- Using Eq. we obtain 

n— 1 

U[ = 2Y,\Pk){Pk\-I, (29) 

fc =0 

where \flk) = RM\flk)- Using Eq. (1^, we obtain 


(30) 


f-|/3A:), if/3fenM = 

^ otherwise, 


where f{k) = 0 if /c G M and f{k) = 1 otherwise. Therefore Um = U[ Uq and using 
Eqs. (29) and (30) it is straightforward to verify that [/( is an orthogonal reflection. ■ 

Next proposition shows that if we mark all vertices inside some polygons of tessellation 
a or fl, the walker is not able to find a marked polygon. 


Proposition 6.2 If the set of marked vertices is a union of polygons of tessellation fl, then 
Um = U. If the set of marked vertices is a union of polygons of tessellation a, then Um is 
similar to U. 
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Proof If the set of marked vertices is a union of polygons of tessellation /3, Eq. (30) reduces 
to 


|«/\ ^ if /5fcnM = 

^ otherwise. 


(31) 


Substituting into Eq. (29), we obtain that U[ = Ui. Therefore Um = U. 

The proof of the second statement is analogous. Define Uq = RmUqRm- Using Eq. ([^, 
we obtain 


m —1 




(32) 


k=0 


where If the set of marked vertices is a union of polygons of tessellation 

a, then using Eq. (Q we obtain 


I ifafcnM = 

“fc/ = S I \ 

Mafc), otherwise. 


(33) 


Therefore Un = Uq and Um = Rm U Rm- 


Useful properties of the QW dynamics are obtained from the spectrum of Um- The 
eigenvalues and eigenvectors can be obtained from the singular values and vectors of the 
discriminant (or Gram-like) matrix D the entries of which are 


Dkk' = 


(34) 


where \oLk) is given by Eq. (Q and |/3^/) is given by Eq. (30). Let and be the 
right and left singular vectors of D respectively and let cos Oj be the corresponding singular 
values obeying = cos0j|/ij) and = cos 0j\vj) and 0 < Oj < 7r/2. Suppose 

that 01 ,..., 9s are the angles in the open interval (0,7r/2), where s is the number of singular 
values of D in the open interval (0,1). Define 


m —1 

k=0 

n—1 

B = ^\l3'k,){k'\. 

k'=0 


(35) 

(36) 


Let A,B be the vector spaces spanned by vectors \ak) and vectors 1/?^/) respectively. We 
reproduce a theorem from [1] adapted to our context. 

Theorem 6.3 The spectrum of Um obeys: 

1. The eigenvalues o/Um with nonzero imaginary part are exactly fgr j = ...,s. 

The corresponding normalized eigenvectors are 


\/2 sin 6 j 


( 37 ) 
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S. An + A'^ nB is the — 1 eigenspace of Um■ A n B^ is spanned by the left singular 
vectors of D with singular value 0. A^ nB is spanned by the right singular vectors of 
D with singular value 0. 

3. AnB + A^ nB-^ is the +1 eigenspace ofUM- AnB is spanned by the left (or right) 
singular vectors of D with singular value 1. 


The above theorem can be used to find most eigenvectors of the evolution operator 
when the singular values and vectors of the discriminant matrix can be explicitly found. 
The eigenvectors that span A'^ H B^ are not described in Theorem 6.3 The staggered 


Fourier transform is an alternative method to find the spectrum of Um, but it depends on 
translation symmetries of the tessellating process. The staggered Fourier transform was 
used in Refs. 


7 Conclusions and Discussions 

In this paper we have formally defined the staggered QW model on graphs. This model 
can be obtained by tessellating the graph using polygons obeying locality restrictions. The 
polygons are cliques and each tessellation covers all vertices. All edges must be inside a 
polygon of some tessellation. An alternative equivalent dehnition is based on the notion of 
orthogonal reflections in the basis associated to the graph, which are unitary and Hermitian 
operators with the following property: The nonzero amplitudes of the eigenvectors associ¬ 
ated with eigenvalue 1 do not have overlap and the sum of those eigenvectors does not have 
zero entries. The evolution operator of staggered QWs must be a product of orthogonal 
reflections that covers all edges of the graph. 

We have showed that any instance of Szegedy’s QW model on a bipartite graph F 
is equivalent to an instance of the staggered QW model on the line graph of F with a 
special kind of tessellations: There is only one vertex in each polygon intersection and 
the tessellation union is a Krausz partition. On the other hand, there are instances of 
the staggered model with two tessellations that cannot be cast into Szegedy’s framework, 
which are the ones with two or more vertices in common in at least one polygon intersection. 
This characterization shows that Szegedy’s model is a restricted version of QWs, which is 
defined only on line graphs of bipartite graphs, a subset of perfect graphs. This fact helped 
researchers to obtain generic results using Szegedy’s model, different from the results about 
the coined version, which is usually found in specific graphs. One of the main results of 
Szegedy’s model is that the quantum hitting time is quadratically smaller compared with 
the classical hitting time on symmetric and ergodic Markov chains. Notice that the classical 
random walk is dehned on a graph F while the quantum walk is defined on the line graph of 
the bipartite graph associated to F. When the expressions for the hitting time on different 
graphs are compared, they have the same structure and it is possible to reach the result. 

There are two main forms of defining searching algorithms in the staggered model: 
(1) Searching d la Szegedy is the method that uses partial tessellations in the generalized 
staggered model. The marked vertices are the ones in the missing polygons. (2) Searching 
d la Falk uses the standard tessellation procedure but interlaces the orthogonal reflections 
with a reflection around the marked vertices. In both cases we can use Szegedy’s spectral 
theorem to calculate the spectrum of the evolution operator when using two tessellations. 
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We have also defined the quantum hitting time on staggered QWs using reflections around 
the marked vertices which generalizes Szegedy’s original definition. We have shown that 
the evolution operator, originally defined as a product of four reflections, can be written as 
a product of two reflections. 

Because all vertices of the two-dimensional lattice belong to the intersection of four max¬ 
imum cliques, it is necessary to employ at least four tessellations. The case with exactly four 
tessellations has one vertex in the polygon intersections but cannot be cast into Szegedy’s 
framework because it employs more than two tessellations. Future works on this issue are 
important to establish the properties of staggered QWs on two-dimensional lattices, the 
analysis of which is missing in literature. 

An interesting subject is the equivalence between the staggered model and the coined 
quantum walk model. This has been established for the line and cycle [Ellis]. A systematic 
analysis on the equivalence of those models on generic graphs is an important issue for future 
works. 
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Appendix 


In this appendix we review the main graph theory concepts [26] we are using in this work 
and we give some examples. 

A simple or undirected graph r(I/, E) is defined by a set V of vertices or nodes and a set 
E of edges so that each edge links two vertices and two vertices are linked by at most one 
edge. Two vertices linked by an edge are called adjacent. The complete graph is a simple 
graph in which every pair of distinct vertices is connected by an edge. A directed graph is 
a graph whose edges have a direction associated with them. 

A subgraph T'{V',E'), where V' C V and E' C E, is an induced subgraph of T{V,E) 
if it has exactly the edges that appear in T over the same vertex set. If two vertices are 
adjacent in T they are also adjacent in the induced subgraph. A clique is a subset of vertices 
of a simple graph such that its induced subgraph is complete. A maximal clique is a clique 
that cannot be extended by including one more adjacent vertex. A maximum clique is a 
clique of maximum possible size. A clique can have one vertex. 

A bipartite graph is a simple graph whose vertex set V is the union of two disjoint sets 
X and X' so that no two vertices in X are adjacent and no two vertices in X' are adjacent. 

A line graph of a simple graph T (called root graph) is another graph T(r) so that each 
vertex of T(r) represents an edge of T and two vertices of T(r) are adjacent if and only if 
their corresponding edges share a common vertex in T. For example, consider the bipartite 
graph r of Fig.j^ The sets of vertices X and X' have labels {ao, ai, 02 , 03 } and {/3o, Pi, 1^2} 
respectively. The edges of F have labels from 1 to 6. Each edge becomes a vertex in the 
line graph T(F) as can be seen in Fig. Now we have to check which vertices of T(F) are 
adjacent. For example, edges 1, 2, and 3 share a common vertex oq of F. They are adjacent 
in T(r) and form a clique with three vertices. This clique has label ao in This goes 

on until we have the partition described in Fig. 



Figure 3: A bipartite graph F and its line graph L{r) with a two-colorable clique partition labeled 
by the vertices of F. 


Szegedy’s model in F uses the following vectors 

|0o) = (|/^o) + |/3i) + 1 ^ 2 )), 

|0i) = |ai)|/3o), \4>2) = \a 2 )\Pi), \h) = \a3)\P2) 
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and 


I 

I'lpo) = (|ao) + |ai)) |/3o), 

|V^i) = ^(l«o> + |a2)) |/3i), 

IV’ 2 ) = ^ (|ao) + las)) \h)- 

The staggered model in T(r) uses the following vectors 

l"o) = (|l) + |2) + |3)) , 

|ai) = |4), \a2) = |5), jaa) = |6) 


and 

|ft) = T(|i) + |4)). 

|A) = T(|2> + |5)), 

|ft) = ^(|3> + |6». 

Notice that there is a bijection between those vectors if we use the following edge-vertex 
correspondence: ao/3o ^ Ij «o/5i ^ 2, ao/32 ^ 3, a;i/3o -H- 4, a2/3i -f-)- 5, 0 : 3/32 o 6. 

If it is given two-colorable tessellations of a line graph with one vertex in each polygon 
intersection, it is possible to find the bipartite root graph. For example, take polygon oq in 
L(r). Selecting the polygons of tessellation /3 that have overlap with oq, we obtain vector 

| 0 o) = (l^o) + 1^1) 1^2)) , 

and the same with respect to polygons ai, a 2 , and 0:3. Now, if we take /3o, polygons ao and 
cci have overlap with /3o, then 

|V’o> = ^(|a:o> + |a:i)) |/3o), 

and the same with polygons /3i and (32- Following this method we re-obtain the vectors 
used by Szegedy’s model in the first place and we re-rebuild the bipartite graph. 

A Krausz partition splits the edges of the graph into cliques in such a way that no 
vertex belongs to more than two of the cliques. For example, the union of the red and 
blue tessellations of graph T(F) in Fig. is a Krausz partition. Notice that each edge is 
in one member of the partition and each vertex is exactly in two members of the partition. 
Krausz [25] proved the following 

Theorem 7.1 A graph F' is a line graph of some graph if and only if F' has a Krausz 
partition. 
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Figure 4: A line graph F' with the Krausz partition and its root graph F. Notice that the partition 
of F' is not two-colorable and F is not bipartite. 


Fig. I shows a graph F' that has a Krausz partition. Then F' is a line graph of some 
graph F, which is depicted in the right side of the figure. Notice again that each edge is in 
one member of the partition and each vertex is exactly in two members of the partition. If 
the Krausz partition is not two-colorable, then root graph is not bipartite. 

In order to define a staggered QW with two tessellations in graph F' of Fig. |^we must 
put two vertices in at least one polygon intersection. Fig. shows that F' is two-tessellable. 
In the simplest case we can take uniform vectors. The vectors of tessellation a (blue) are 

l“o) = 2 (|l) + |2) + |3) + |4» , 

|ai) = |5) 

and the vectors of tessellation f3 (red) are 

|A) = A{|i> + |2)), 

bi) ^ {|3) + |4) + |5)) . 

It is interesting to analyze staggered QWs on graphs that are not line graphs. Those 
graphs do not have root graphs and any comparison with classical random walks must be 



Figure 5: Tessellations of graph F'. 
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Figure 6: Nine forbidden Beineke induced subgraphs. 


performed on the same graph. Krausz partitions play no role in this case. There is an 
alternative form to characterize line graphs. Beineke [23] proved the following 

Theorem 7.2 Let T' be a graph. There exists a graph T such that T' is the line graph ofT 
if and only if T' contains no graph of Fig. as an induced subgraph. 

Graph Ti of Fig.[^is called claw. A trivial consequence of Beineke’s theorem is the following 


Corollary 7.3 A line graph is claw-free. 

By inspection, we verify that graphs r 2 , Ta, r 4 , and Tg can be tessellated with two 
tessellations if we put two vertices in common in a polygon intersection while Ti, Tg, r 7 , 
Tg, and Tg cannot be tessellated with two tessellations. Tg cannot be tessellated with 
two tessellations because of parity violation and the remaining graphs have three or more 
maximal cliques with a common vertex. To define a staggered QW on a graph that is not 
a line graph of a bipartite graph, either we have to use a tessellation that has at least two 
vertices in common in a polygon intersection or we have to use more than two tessellations. 

As a last example, consider the two-tessellable graph of Fig. which is not a line 
graph as one can check by finding one of the Beineke induced subgraphs. The vectors of 
tessellation a (blue) are 

l“o) = 2 (|0) + |l) + |2) + |3)) , 

K) = (W +15» 

and the vectors of tessellation fd (red) are 

|A) = T(|o> + |i)). 
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Figure 7: A two-tessellable graph that is not a line graph. The tessellations have two vertices in 
common in at least one polygon intersection. 


|/^i) - 2 (1^) + 1^) 1^) 1^)) • 


This choice generates the unitary operator 



1 -1 

-1 1 

0 0 

0 0 

1 1 

1 1 


1 1 

1 1 

1 -1 

-1 1 

0 0 

0 0 


0 0 

0 0 

1 1 

1 1 

1 -1 

-1 1 


which is obtained from Eqs. Q to ([^. U satisfies = I. Therefore the evolution is 
periodic. To obtain a more general evolution operator we can choose |/3o) as a non-trivial 
point in the Bloch sphere by taking 


|/3o) = cos^ |0) + e *^sin^ |l), 


where 0 < 9 < tt and 0 < (/? < 2tt. 

Two-tessellable staggered QWs are interesting because we can use Szegedy’s spectral 
theorem to find the spectrum of the evolution operator. All line graphs of bipartite graphs 
are two-tessellable and generate staggered QWs equivalent to Szegedy’s QWs. Not all line 
graphs of non-bipartite graphs are two-tessellable and the same is true for graphs that 
are not line graphs of any graph. It would be interesting to characterize the class of two- 
tessellable graphs. 


21 












